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(§4) Precede de trartement d'fmage de haut niveau, dispositif pour sa mise en oauvre et son application a la de- 
tection d'objets ponctuels. 

fen Le procede consiste dans une premiere etape (1), a 
calculer une fonction de coOt (C) a partir d'erreurs issues 
de la comparaison entre une detection d'objets dans un lot 
d'images determine resultant du trartement de haut niveau, 
et un resultat d'une fonction de detection predeterminee ( 
5 J associee au meme tot d'images, appelee verite-terrain, ; 
et consiste, dans une deuxieme etape (2); pour minimiser 
iterativement par approches successives la fonction de 
coOt (C), a initialiser un parametre d'optimisation (T) a une 
valeur nnie determinee (T 0 ) survant une repartition determi- 
nes du coOt minimal (C^^J, le processus iteratjf s'arretant 
quand le minimum globaiae la fonction de coOt C(T) f cor- 
respondant au reglage parametrique optimal, est atteiht 

Applications: analyse d'image et, notamment, detection 
d'objets ponctuels. 
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Proced<§ de traitement d'image de haut niveau, dispositif pour sa mise 
en oeuvre et son application h la detection d'objets ponctuels. 

5 

L'invention concerne un precede de traitement d'images de haut 
niveau, un dispositif pour sa mise en oeuvre et son application a la detection 
d'objets ponctuels. 

Uinvention se rapporte au domaine de I'analyse d'images. 

10 Dans les traitements d'images, un traitement d'image de haut 

niveau permet d'obtenir une collection d'objets symboliques a partir d'objets 
detectes dans une image pour faire de la detection, de la reconnaissance de 
forme, etc... Une chaTne de traitement de haut niveau doit etre reglee eh 
ajustant un jeu de parametres qui conditionnent ses performances. 

15 Pour obtenir un reglage parametrique, une premiere methode 

connue consiste a realiser manueflement des essais successifs jusqu'a 
I'obtention de resultats juges satisfaisants. Une deuxiemfc methode consiste 
a mettre en evidence qu'une erreur de traitement est une fonction analytique 
des parametres. Elle consiste ensuite a minimiser I'erreur par une methode 

20 de descente du gradient. Cette deuxieme methode est notamment utilisee 
pour calculer les poids des connexions d'un reseau de neurones. 

Cependant la premiere methode dite m§thode manuelle, 
subjective, est rapidement inapplicable des que la nombre de parametres 
augmente et la deuxieme methode presuppose I'existence d'un lien 

25 analytique avec I'erreur. 

Le but de l'invention est de pallier les inconvenients precites. 
A cet effet, l'invention a pour objet un precede de reglage des 
parametres d'une chaTne de traitement d'image de haut niveau utilisant un 
processus iteratif d'optimisation, caracterise en ce qu'il consiste, dans une 

30 premiere etape, a calculer une fonction de cout a partir d'erreurs issues de 
la comparaison entre une detection. d'objets dans un lot d'images determine 
resultant du traitement de haut niveau, et un resultat d'une fonction de 
detection predetermines associee au meme lot d'images, appelee verite- 
terrain, et en ce qu'il consiste dans une deuxieme etape t pour minimiser 

35 iterativement par approches successives la fonction de cout, a initialiser un 
parametre d'optimisation a une valeur finie deterrninee suivant une 
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repartition determinee du cout minimal, a appliquer des perburbations 
aleatoires, pour chaque valeur du parametre d'optimisation dont la variation 
obeit a une fonction recursive decroissante determinee, sur les parametres 
de traitement pendant un nombre determine ^iterations, et si les 
perturbations sont acceptees selon un critere determine, a remplacer les 
parametres courants de traitement par les parametres apres perturbations 
pour le reglage des parametres, sinon a conserver les parametres courants 
de traitement, et a tester si la valeur minimale de la fonction de cout repond 
a tin critere de convergence determine, si le test est negatif, a recalculer la 
valeur du parametre d'optimisation en suivant la fonction recursive 
decroissante determinee, et si le test est positif, a arreter I'etape de 
minimisation, le minimum global de la fonction de cout correspondant au 
reglage parametrique optimal etant atteint 

Uinvention a pour avantage qu'elle permet d'une part le reglage 
automatique par minimisation d'un coQt de detection entre le resultat du 
traitement et un resultat attendu appele verite-terrain, en ne faisant aucune 
hypothese sur le fonctionnement interne du traitement et permet tfautre part, 
la convergence du cout vers un minimum global, meme en presence de 
minima locaux. De plus, le processus it6ratif de minimisation du cout, 
mesurable en tout point d'un espace d'etat determine, ne requiert aucune 
definition analytique du cout en fonction de I'etat 

D'autres caracteristiques et avantages de la presente invention 
apparaitront plus clairement a la lecture de la description qui suit, faite en 
regard des figures annexees qui represented : 

- la figure 1, un organigramme des etapes du precede selon 

Tinvention, 

- les figures 2a, 2b et 2c, des representations symboliques 
respectives a une bonne detection, une fausse alarme et £ une non- 
detection, 

- la figure 3, un organigramme de la phase ^initialisation du 
parametre d'optimisation T f 

- la figure 4, un organigramme de I'etape de minimisation du cout 
du precede selon Tinvention, 

- la figure 5, un schema de principe d'un mode de realisation d'un 

' dispositif pour la mise en oeuvre du procede selon Tinvention, et 
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- la figure 6 ( un exemple d'application numerique presents sous la 
forme d'un tableau. 

Le procede selon I'invention est illustre par I'organigramme de la 

figure 1 . 

5 II consiste dans une premiere etape 1 a calculer une fonction de 

cout d'erreurs de detection C, en comparant le resultat d'un traitement de 
haut niveau sur un lot d'images representatif d'un objectif determine avec un 
resultat attendu issu d'une fonction de detection pr6determinee, appel6e 
"verite-terrain", sur le meme objectif. La fonction de cout C tient compte de la 

10 nature des erreurs de detection ainsi que des erreurs sur les mesures, en 
les ponderant pour tenir compte des besoins operationnels. 

Le procede se poursuit dans une deuxieme etape 2 par la 
minimisation du cout C(Ti) en fonction d'un parametre d'optimisation 71 pour 
le reglage, par approches successives, des parametres P du traitement. 

15 Comme decrit precedemment, le traitement d'images de haut 

niveau permet d'obtenir une collection d'objets symboliques a partir d'objets 
detectes dans une image. Les objets symboliques, pour etre identifies et 
caracterises, sont munis d'une representation sous forme de vecteur de 
mesure, de liste, de chatne ou de graphe. II est possible de quantifier les 

20 similitudes entre ces objets en calculant une "distance" entre ces 
representations. 

Cette distance est metrique ou binaire dans le cas de vecteurs. 
Elle correspond a un cout de mise en correspondence entre chaines, appele 
egalement distance d'edition qui est calculee a partir de I'algorithme de 
25 Wagner et Fisher, ou correspond encore a un cout de mise en 
correspondance entre graphes. 

Dans la presente invention les objets detectes sont representes 
par des vecteurs de mesure de dimension donnee m. Les mesures sont des 
nombres reels et Pespace des mesures est identifie a R m . La fonction de 
30 detection 5 utilisee par le traitement est definie de la maniere suivante : 
8:i^P(R m ) 
(1) I->6(I)={ M ,M..»} 
ou i est I'ensemble de depart representant un lot d'images et P(R m ) 
['ensemble d'arrivee representant un ensemble de parametres P definissant 
I'espace des mesures R m . La fonction de detection 5 definie sur i fait 
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correspondre a une image I un nombre de vecteurs de mesure 
5(l)={ni,H2-}- 

La verite-terrain est definie comme une fonction de detection 
particuliere 8 V connue a priori. 
5 Elle correspond, par exemple, a une representation d'avions 

detectes sur une piste delimitant un objectif determine, ou zone, a partir de 
laqueile sont determines I'abscisse et I'ordonnee du centre d'inertie de la 
zone. 

Chaque image du lot d'images est traitee avec un reglage 
10 parametrique courant fonction d'un contexte determine dans lequel 
fonctionne la chaThe detraitement appele egalement etat du systeme. 

Un objet detecte est done caracterise par un vecteur de mesure \i 
et un objet de la verite-terrain par un vecteur v 

La premiere etape 1 du procede consiste dans une premiere 
15 phase, pour calculer la fonction de cout C, a associer de maniere unique un 
objet detecte avec un objet de la verite-terrain caracterises respectivement 
par leur vecteur uetv. Cette association se fait en respectant deux criteres 
successifs. Dans un premier temps les vecteurs de mesure nicaracterisant 
les objets detectes et les vecteurs vj caracterisant les objets de la verite- 
20 terrain doivent verifier le critere suivant : la valeur absolue de la difference 
entre la i-eme coordonnee du vecteur de mesure n et la i-eme coordonnee 
du vecteur verite-terrain v, soit I'ecart entre les mesures, doit etre 
strictement inferieure a urie valeur oq determinee traduisant une incertitude 
maximale sur les mesures. Cette valeur permet de traduire les differences 
25 d'ordre de grandeur ou de sensibilite entre les mesures. Cependant, un objet 
detecte peut etre associe a plusieurs objets de la verite-terrain. 

Une telle ambigutte est .levee dans un deuxieme temps en 
affectant preferentiellement un objet detecte a I'objet de la verite-terrain le 
plus proche au sens de la distance de Hamming ponderee donnee par la 
30 formule suivante definissant le deuxieme critere dissociation : 

m-1 

(2) d a (u,v)= Elwc-vkKak 
k=0 

Le nombre d'associations, c'est-a-dire le nombre de couples 

(Hi ,v;) de vecteurs verifiant les deux criteres precedents, definit un nombre 

de bonnes detections N bd . 
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Les vecteurs "celibataires", c'est-a-dire non associes suivant les 
deux criteres precedents, sont cpnsideres comma des erreurs de detection. 

Les erreurs de detection sont de deux natures differentes definies 
respectivement par un nombre de fausses alarmes et de non-d6tections 
5 N nd . 

Dans un exemple d'un lot d'images comportant un nombre N(k) 

determine d'images et des objets de la verite-terrain associes, & chaque 
image k est denombre un nombre N v (k) determine d'objets de la verite- 
terrain. A Tissue du traitement de detection sur I'image k, param6tr6 par un 
10 jeu determine de parametres, un nombre determine N<j(k) d'objets detectes 

est obtenu. Dans une deuxieme phase le precede recherche le nombre de 
bonnes detections N^k) pour I'image k. Le nombre de fausses alarmes 

N fe (k) et le nombre de non-detections N nd (k) sont done calcules par les 
expressions suivantes : 
15 (3) N nd (k) = N v (k)-N bd (k) 

(4) N fe (k)=N d (k)-N bd (k) 

Les figures 2a, 2b et 2c illustrent des representations 
symboliques de "bonnes" et de "mauvaises" detections, dans un espace de 
mesure de dimension egale a deux ; les mauvaises detections etant 
20 traduites par des non-detections et/ou des fausses alarmes. Sur les trois 
figures, des objets detectes 4, 4^ , 42 sont representes par un petit cercle, un 
objet associe de la verite-terrain 5 par un petit rectangle et une zone 
d'incertitude de mesure 6 est representee par un rectangle en trait 
discontinu. 

25 Pour qu'il y ait une "bonne detection", il faut que I'objet detecte 4, 

4i , 42 et I'objet de la verite-terrain 5 associe se retrouvent au plus proche 
Tun de Tautre au sens de la distance de Hamming ponderee a Tinterieur de 
la zone d'incertitude de mesure 6. 

La figure 2a represente une "bonne detection". Un objet unique 

30 detecte 4 se trouve a proximite de I'objet de la verite terrain 5 dans la zone 
^incertitude de mesure 6. Les deux criteres dissociation sont verifies pour 
I'objet detecte 4. 

La figure 2b represente une "bonne detection" similaire a celle de 

la figure 2a et deux fausses alarmes correspondant respectivement a un 
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eloigne de I'objet de la verite-terrain 5, et a un deuxieme objet detects 42 
situe en dehors de la zone d'incertitude de mesure 6. Le deuxieme critere 
donne par la formule (2) n'est pas verifie. 

La figure 2c represent? une "non-detection". Aucun objet detecte 
5 4 ne se trouve dans la zone d'incertitude de mesure 6. Aucun critere n'est 
verifie. 

Selon une troisieme phase, le procede calcul pour un lot 
determine d'images un cout d'erreur de detection C a partir d'un 
pourcentage de non-detection P nd , d'un pourcentage de fausse-alarme P fa 
10 et d'une erreur relative moyenne D bd sur la mesure des objets correctement 
detectes. 

Les pourcentages de non-detection P nd et de fausses alarmes 
P fe sont definis respectivement par les formules suivantes : 
nd Ef_ 0 N nd (k) 

Zk_ 0 N v (k) 



15 



20 



25 



(6) P fe = 



Z^oN & (k) 
Ek=0 N d( k ) 



L'erreur relative D bd correspond a la distance moyenne entre les bonnes 
detections et les objets associes de la verite-terrain. Cette distance est 
calculee par la formule suivante : 

bd 1 £ ^ = o2^» v associes - correctement d a (n, v) 

(7>D = ™ E^ =0 Nv(k) 
Le cout d'erreur de detection C est considere comme une combinaison 
lineaire a coefficients positifs et de somme egale a 1, du pourcentage de 
fausses alarmes P fe , du pourcentage de non-detection P na et de l'erreur 
relative moyenne sur les mesures D bd . II est donne par la formule suivante : 

(8) C = M-P fa +k2-P nd+ *3- Dbd 

Les coefficients de ponderation Xi.XgetA* sont determines de maniere a 
traduire I'importance relative accordee a chaque type d'erreur en fonction 

des besoins operationnels et verifient la relation suivante : 

(9) X 1+ *2+*3= 1 



7 
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Le procede consiste ensuite dans une deuxieme etape 2 a 
minimiser iterativement la fonction de cout C precedemment calculee dans la 
formule (8) en utilisant une melhode iterative d'optimisation appelee 
methode de recuit simulS adaptee au reglage parametrique d'une chaTne de 
5 traitement d'images de haut niveau. 

Une description de cette methode peut etre trouvee dans un 
article intitule : "Simulated annealing : Theory an application" de P.J.M Van 
Laar Hoven et E.H.L Aaarts et publie par D. Reidee Publishing Company 
1987. 

t 

10 La methode de recuit simule tire son nom de la thermodynamique 

statistique: lorsque la temperature d'un ensemble de particules en 
interaction dans un systeme est abaisse de facon suffisamment lente, celui- 
ci atteint un etat cfenergie minimale. 

A chaque etat E p du systeme est associe un cout C p assimilable 

15 a une energie. La methode modifie ateatoirement I'etat E D en I'etat E n de 
coflt C q . Soft AC pq =C q -C p . P q 

Si ACpq ^ 0, la modification ou perturbation est acceptee. Si le 
cout augmente, la perturbation peut etre acceptee avec la probabilite 
exp(-AC pq /t) ou T est la "temperature du systeme". L'interet d'accepter les 

20 augmentations de cout est de permettre au systeme de converger vers un 
minimum global alors que la methode du Gradient Stochastique, n'acceptant 
que des diminutions de cout, peut faire converger le systeme vers un 
minimum local. 

Les divers travaux effectues sur I'application du recuit simule 
25 exposent une strategie de refroidissement optimal dans le cas ou la 
repartition du cout est gaussienne autour du cout moyen. Or, le cout lors de 
calibration d'algorithme n'a pas necessairement une repartition gaussienne. 

Le procede selon Pinvention utilise une "strategie de 
refroidissement" adaptee, en modifiant I'initialisation du parametre 
30 d'optimisation. Une caracteristique du cout de detection C est qu'il reste 
eleve pour une majorite de reglages parametriques sans beaucoup varier. 
Par contre, la partie de I'espace parametrique dans laquelle le cout est 
proche du cout optimal est relativement petite. Loin de cette region, il est 
difficile de la rejoindre, bien que statistiquement cela est toujours possible. 
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Pour minimiser le temps de calcul, le procede recherche d'abord 
cette region interessante on, en fait, de bonnes detections sont enregistrees, 
puis s'efforce d'y rester, tout en continuant la recherche du minimum. 

Pour faciliter la comprehension et conserver une certaine 
5 cohesion dans le raisonnement de la methods du recuit simule adaptee au 
reglage des parametres d'une chaine de traitement de haut niveau, la 
description du procede selon I'invention conserve I'analogie avec le terme 
"temperature T" qui devra etre considere comme un parametre 
d'optimisation determine. 
10 Le procede seton I'invention consiste au debut de la deuxieme 

etape 2 dans une premiere phase a initialiser le parametre d'optimisation T. 

La temperature initiate doit etre suffisamment elevee pour 
"oublier" une configuration initiale de la chaine de traitement, par example 
un lot d'image precedent definissant un objectif de nature differente. Cette 
15 temperature peut etre fixee arbitrairement ou fixee en fonction de la 
repartition du cout C a cette temperature. La figure 3 illustre la phase 
^initialisation du parametre T. 

Celled a lieu sur un npmbre determine A d'iterations visant a 
rechercher le cout le plus faible C m jn-init Pour chaque iteration L. une 
20 etape 7 calcule le cout Cl. suivant la troisieme phase de la premiere etape 1 
puis, dans une etape 8, des perturbations sur les parametres Pl de 
traitement sont appliquees a la chaine de traitement. Dans cette phase 
d' initialisation toutes les perturbations sont acceptees. Une etape 9 calcule 
I'effet des perturbations sur le cout Cl et une etape 10 ajuste les parametres 
25 Pl du traitement en fonction du resultat de I'etape precedente 9. Un test 11 
verifie si le nombre determine d'iterations A est atteint Si le resultat du test 
11 est negatif, un compteur d'iterations, L, 12, est increments et le 
processus continue. Si le resultat du test 1 1 est positif, le meilleur cout c'est- 
a-dire le cout minimal C m j n _j n jt et les parametres associes de traitement Pl 
30 sont conserves dans une etape 13. Le cout C m jn-init 61 le i eu de 
parametres correspondant Pl devierinent les valeurs initiates pour la 
premiere iteration a une temperature finie To qui est determinee de la 
maniere suivante : 

un deuxieme test 14 compare le cout minimal C m j n _j n it avec le 
35 coefficient X2 ponderant les erreurs de non-detection. Si 0^.^ <k 2 , il 
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existe necessairement de bonnes detections et la temperature finie Tq est 
fixee de maniere a minimiser les chances de "remonter" vers une zone de 
cout eleve ou il n'y a pas de bonnes detections. Une probabilite determinee 
P-l d'accepter une augmentation du cout C m j n .j n jt est fixee pour Pf <0,5. La 
5 temperature finie Tq est alors determinee par la formule suivante dans une 
etape 15: 

(10) To=-ln(P 1 ).C rnin -i nit 

Sinon, une probabilite cFaccepter un accroissement de cout de X1+A.3 de 
maniere a continuer la recherche de la zone de faible cout sans raster 
10 bloque par une barriere, est fixee pour P2<0,5. Le cout maximal 6tant <§gal a 

1 1 AC max = 1 ~ ^2 = *4 + ^3- 

Dans ce cas la temperature finie Tq est determinee par la formule 
suivante dans une etape 16 : 

(11) 1]}=-In(P2).(^1+^3) 

15 La figure 4 illustre une deuxieme phase de la deuxi&me Stape 2 

du procede selon I'invention. 

La temperature finie Tq venant dletre determinee dans la 

premiere phase de la deuxieme etape 2 du procede, (a temperature T 

decroTt a chaque iteration i, suivant une loi de decroissance determinee, en 
20 fonction de la variation de la fonction de coQt C. La temperature pourrait 

varier selon la loi geometrique suivante : 

(12) 1+1 1 
ayec 0<a <1 

Cependant, le cout ne varie pas uniformement lorsque la temperature 
decroTt. Des phenomenes de transition de phases peuvent etre observes 

25 lorsque le cout moyen varie fortement. II est opportun de diminuer lentement 
la temperature lors de ces transitions de phases afin de reduire les risques 
de piegeage dans un minimum local. Par centre, lorsque le cout moyen varie 
peu ( le syst£me est dans une phase stable et il convient d'en profiter pour 
faire varier fortement la temperature. 

30 La loi de decroissance de la temperature utilisee par le procede 

selon T invention tenant compte de ce comportement est donnee par la 
formule suivante : 

(13) T l + 1 =T Lex p[-d.T I /a(T 1 )] 
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Ou d est une constante determinee et ou o(Tj) correspond a I'ecart type sur 
le cout C a. la temperature Ti 

Un compteur d'iterations i , est initialise a zero. II est ensuite 
increments a chaque iteration i sur la temperature Tj. L'etape 17 calcule a 
5 chaque iteration i une temperature T,+i a partir de Iteration precedente 
^ suivant la loi (13) de decroissance definie precedemment. 

L'etape 1 8 consiste a perturber I'etat E p du systeme associe a un 
cout C p pour obtenir un nouvel etat Eq associe a un nouveau cout C q . A 
chaque iteration i , les parametres p(Ti) sont tous modifies aleatoirement 
10 puis un a un, ce qui reduit artificiellement la dimension de I'espace d'etats 
done le nombre d'iterations i necessaires. La nouvelle valeur d'un 
parametre est choisie dans un voisinage de la valeur courante avec une loi 
uniforme et un ecart maximal qui est ajuste en fonction de revolution du cout 
C(Tj). La perturbation des parametres P(Ti) est donnee par la formule 

15 suivante : 

hav ad (t V ap fT \ <C(Ti-l)>-KC(Ti)> 
(14) AP max (T i+ i j = AP max (Tij. 2<C ( T .) > 

ou < C(Ti_i) > et < C(Ti ) > repr§sentent respectivement la valeur moyenne du 
cout a la temperature de I'iteration precedente et la valeur moyenne du cout 
a la temperature de Iteration courante. 
20 Ainsi les parametres sont d'autant moins perturbes que 

APmaxCTi+l) converge vers une solution. La valeur initiate de la perturbation 

maximale des parametres est fixee par I'utilisateur. 

L'etape 19 consiste a calculer une variation du cout AC pq definie 

comme la difference entre le cout C q associe a I'etat Eq, correspondant au 
25 nouvel etat du systeme apres perturbation, et le cout C p associe a I'etat E p 
correspondant a I'etat du systeme ayant perturbations. Un premier test 20 
teste si la variation du cout AC pq est inferieure ou egale a zero, e'est-a-dire 

si le nouveau cout C p est inferieur ou egal a I'ancien cout C p traduisant une 
diminution du cout. Si le test 20 est positif alors la perturbation est acceptee 
30 et un nouveau jeu de parametres P(Ti)est genere pour le reglage des 

parametres de traitement represente par l'etape 21. Si le test 20 est negatif 
e'est-a-dire si la variation AC pq est superieure a zero, un deuxieme test 22 

permet d'accepter la perturbation si le resultat du test 22 repond a la 

probabilite suivante : 
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(15) exp(-AC pq /Ti) > random[0,l] 

c'est-a-dire que la fonction exponentielle du rapport de la variation negative 

du cout avec une temperature determinee du systeme doit etre sup6rieure & 

un tirage aleatoire compris entre 0 et 1 et les parametres de reglage sont 
5 modifies de la valeur de la variation AC pq a I'etape 21. Si le resultat du test 

22 est negatif le nouveau jeu de parametres est rejete et une nouvelle 

perturbation sur les parametres de traitement precedents est g§ner§e. Le 

processus est reitere tant que le systeme n'a pas atteint Tequilibre 

thermodynamique approche" terme utilise par la methode de recuit simule. 

10 Pour le reglage de parametres, il suffit de choisir un nombre determine 
Lj(Ti) d'iterations j pourchaque temperature Tj. 

Le nombre d'iterations effectuees a une temperature 1} est donne 
par la formule suivante : 

(16) Lj(T 1 ) = A.«p(b.o(T l -i)/T 1 ) 

15 A est un coefficient de probability choisi egal au nombre d'iterations de la 
phase d'initialisation a la temperature finie Tq. Lj est limite a 5A pour limiter 

les temps de calcuL b est une constante determinee par I'utilisateur et 
o(Ti-i) est I'ecart type du cout a la temperature de I'iteration precedente 
i-1. 

20 Un troisieme test 23 teste si le nombre d'iterations Lj(li) est 

atteint. Tant que le nombre maximal d'iterations Lj(Tj) n'est pas atteint, le 

compteur d'iterations j s'incremente et les parametres de traitement sont ou 
ne sont pas reajustes pour une meme temperature T\ determinee. Une fois 
le nombre d'iteration maximal Lj(T x ) atteint, un quatrieme test 24 permet de 

25 verifier si le systeme est "suffisamment refroidi" au sens de la methode du 
recuit simule. 

Pour le reglage des parametres, le processus iteratif de 
minimisation s'arrete quand un critere d'arret d§termine est verifie. 

Le critere d'arret du processus est determine par la formule 

30 suivante : 

(17) <C(T i )>-C min ^ 

< Cmin-init > "^min 
II permet d'interrompre le processus lorsque le cout se stabilise a un faible 

niveau. Dans la formule (17) < C(Tj) >et< C m [ n _j n ; t > correspondent 

respectivernent a !a valeur moyenne de la fonction de cout abtenue a la 
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temperature Ti et a la valeur moyenne du cout minimal a I'initialisation, 
correspond a la plus petite valeur obtenue pour la fonction de cout depuis le 
debut du processus d'optimisation. s correspond a un parametre de 
convergence determine par I'utilisateur. Si le critere d'arret n'est pas verifie, 

5 le compteur ^iterations i est increments et une nouvelle temperature 
est calculee suivant la formule (13) jusqu'a ce que le cout minimal converge 
vers le parametre de convergence e. De plus, le procede sauvegarde 
systematiquement le meilleur reglage parametrique et son cout associe pour 
former le resultat de la calibration du processus d'optimisation. 

10 Le minimum global est alors atteint et les parametres sont 

inchanges jusqu'a ce qu'un nouveau lot d'images soit traite par la chaTne de 
traitement de haut niveau. 

La figure 5 illustre un schema de principe d'un dispositif de 
traitement tfimage de haut niveau utilise pour la mise en oeuvre du procede - 

15 selon I'invention. 

II comporte des moyens de traitement 25 recevant sur une 
premiere entree un nombre determine d'images formant un lot d'images. Us 
recoivent sur une deuxieme entree un ensemble de parametres P decrivant 
I'etat du systeme pour chaque iteration. La sortie des moyens de traitement 

20 25 est couplee a une premiere entree de moyens de calcul 26 du cout C qui 
recoivent sur une deuxieme entree des informations issues de la verite- 
terrain. La sortie des moyens de calcul 26, delivrant un cout d'erreur de 
detection pour chaque iteration, est couplee a une premiere entree de 
moyens de minimisation du cout C, 27, realisant les etapes de la methode 

25 de recuit simule adaptee pour le reglage des parametres P du traitement. 

Les moyens de minimisation du cout 27 recoivent sur une 
deuxieme entree les perturbations aleatoires des parametres de traitement 
P. La sortie des moyens de minimisation du cout C, 27, delivre une valeur de 
variation des parametres P en fonction des perturbations aleatoires. Elle est 

30 couplee a la deuxieme entree des moyens de traitement d'images 25. 

Le procede selon I'invention peut etre applique a la detection 

d'objets ponctuels. 

Dans le cas c^une chaine de detection d'objets ponctuels, les 

mesures caracteristiques des objets extraits sont les coordonnees de la 

35 position de I'objet, abscisse et ordonnee dans I'image. L'objet etant pon'ctuel, 
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I'erreur de mesure autorisee est du meme ordre de grandeur que la 
resolution spatiale du capteur, qui correspond a un pixel dans I'image. Seuls 
les fausses alarmes et les non-detections sont representatives. Les 
coefficients de pohderations relatifs au type d'erreurs de detection 

5 repondent alors a une nouvelle relation donnee par la formule suivante : 
(18) (m,X2.*3) = (M-M) 
Pour cette application, les moyens de traitement 25 de la chaine de 
detection d'objet ponctuel comportent des moyens de filtrage pour la 
reduction du bruit, des moyens de seuillage de la luminance de chaque pixel 

10 de ' r image, des moyens de recherche des composantes connexes, des 
moyens de calcul de la taille des objets detectes dans I'image correspondant 
a une aire mesuree en nombre de pixels, et des moyens de selection des 
objets dont la taille est comprise entre deux, valeurs limites determinees. Les 
objets detectes sont dans des regions de I'image de taille tres faible de 

15 I'ordre de quelques pixels, dont la position est obtenue en recherchant le 
centre d'inertie de ces regions. Trois parametres sont alors regies : le seuil 
de detection, la taille minimale et la taille maximale des objets. 

A titre d'application numerique, la figure 6 illustre par un tableau, 
un exemple de valeurs correspondant aux differentes variables accessibles v 

20 par I'utilisateur : 

X represente le coefficient de ponderation des erreurs de 
detection, Pi une probabilite d'accepter une augmentation du cout minimal 
c min> p 2 une probabilite d'accepter une augmentation de X,e le parametre 
de convergence et d et b sont des constantes ponderant I'ecart type du cout 

25 C. 



i4 
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REVINDICATIONS 

1 . Prpcede de reglage des parametres d'une chatne de traitement 
d'image de haut niveau utilisant un processus iteratif d'optimisation, 

5 caracterise en ce qu'il consiste, dans une premiere etape (1), a calculer une 
fonction de cout (C) a partir d'erreurs issues de la comparaison entre une 
detection cfobjets dans un tot d'images determine resultant du traitement de 
haut niveau, et un resultat d f une fonction de detection predeterminee (8 V ) 

associee au meme lot d'images, appelee verite-terrain, et en ce qu'il 
10 consiste dans une deuxieme etapes (2), pour minimiser iterativement par 
approches successives la fonction de cout (C), a initialiser un parametre 
d'optimisation (T) a une valeur finie determinee (To) suivant une repartition 
determinee du coQt minimal (Gmjn-init), a appliquer des perburbations 

aleatoires, pour chaque valeur du param&tre d'optimisation (Tj) dont la 
15 variation obeit a una fohction recursive decroissante determinee, sur les 
param&tres de ^ traitement! (Rj) pendant un nombre determine cfit§rations (Lj), 
et si les perturbations sent acceptees selon un crit§re determine, a 
remplacer les parametres courants de traitement par les parametres apres 
perturbations pouh le reglage des parametres, sinon d conserver les 
20 parametres courants de traitement, et a tester (24) si la valeur minimale 
( c min) de la fonction de pout (Cxi) repond a un crit§re de convergence 
determine, si le test (24) est negatif, a recalculer la valeur du parametre 
d'optimisation (Tj) en suivant la fonction recursive decroissante determinee, 
et si le test (24) est positif, a arreter retape de minimisation (2), le minimum 
25 global de la fonction de cout (Cxi) correspondant au reglage parametrique 

optimal etant atteint 

2. Pr6ced6 selon la revendication 1, caracterise en ce qu'il 
consiste, pour le calcul de Isi fonction de cout (C) r e assbcier, suivant deux 
criteres d'associations determines, un vecteur de mesure repr6sentatif d'un 

30 objet detecte par le traitement avec un vecteur de mesure unique 
representatif d'un objet de la v6rite-terrain ( et caracterise en ce que la 
fonction de cout (C) est une combinaison lineaire, a coefficients (Xj, X2* ^3) 

positifs et de somme egale a Un, des taux d'erreurs de detection (P fa ,P nd ) 

defintes a partir des objets detectes (Nd(k)), dans chaque image (k), des 

35 objets detectes verifiant les criteres d'association (N db (k)) et des objets de 
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la verite-terrain (N v (k)), et de I'erreur relative moyenne (D bd ) sur la 
mesure des objets verifiant les criteres d'association. 

3. Proced6 selon la revendication 2, caract6ris6 en ce que le 
premier critere d'association consiste a selectionner les vecteurs mesures 

5 {it et v), representatifs des objets detectes et des objets de la. verite-terrain, 
dont Pecart de leurs coordonnees respectives et vi)est inf6rieur a une 
valeur determinee traduisant une incertitude maximale (cq) sur les mesures, 
et en ce que le deuxieme critere d'association consiste a affecter 
preferentiellement, parmi les vecteurs selectionnes par le premier critere, un 

10 objet detecte a un objet de la verite-terrain le plus proche au sens de la 
distance de Hamming ponderee. 

4. Procede selon la revendication 2, caracterise en ce que les 
taux d'erreurs de detection pour chaque image (k) sont la somme d'un taux 
de fausses-alarmes (P fe ) correspondant au rapport de la difference entre le 

15 nombre d'objets detectes (Nd(k)) dans Pimage (k) et le nombre d'objets 

detectes verifiant les criteres d'association (N bd (k)) et le nombre cfobjets 

(Nd(k)) detectes dans I'image (k), et e d'un taux de non-detection (P nd ) 

correspondant au rapport de la difference entre le nombre d'objets de la 
verite-terrain (N v (k)j et le nombre d'objets detect6s verifiant les criteres 

20 d'association (N bd (k)) et du nombre d'objets de la verite-terrain (N v (k)), et 
en ce que I'erreur relative moyenne (D bd ) correspond ± la distance 
moyenne entre les objets d§tectes verifiant les criteres d'association et les 
objets de la verite-terrain respectifs. 

5. Procede selon les revendications 1 a 4, caracterise en ce que 
25 Initialisation du parametre d'optimisation a une valeur finie determinee (To) 

consiste a rechercher par un tirage al6atoire un jeu de parametres (Pj) 
permettant d'pbtenir au moins une bonne d6tection et choisi pour borner 
superieurement la probabilite de ne plus avoir de bonnes detections. 

6. Procede selon la revendication 5, caracterise en ce qu'il 
30 consiste a effectuer un nombre determine (A) d'iterations (L) et a calculer (7) 

un cout minimal (C m j n .j n it) ( les parametres (P|J de traitement relatifs a ce 
cout (C m j n .j n jt) initialisant les paramdtres de reglage, et S comparer ce cout 
(C m j n .j n j t ), par un test (14), a un coefficient (A.2) de ponderation 

correspondant aux erreurs de non-detection et si le cout minimal (C m j n _j n j() 



est inferieur a ce coefficient ; (^), a calculer la valeur finie du parametre 
d'optimisation (Tq)}en for^i^^ valeur de probabilite (Pi) et 

du cout minimal XC^^^k^piie^^>nriinimai (C m ir^hjt) est superieur au 
coefficient (A*), a calcu^^leMr^hie du parametre optimisation (Tn) en 
5 fonction d'une deuxieme valeur. de probabilite (P2) et des coefficients de 
ponderation relatifs auxfausses alarmes (^) et aux erreurs sur les mesures 

7 Process sei^;le^ 1 a 6; caracterise en ce que 

les perturbations aleah|res des parametres de traitement (Pr) sbnt 
10 obtenues en calojlant, suivant une loi uniforme determinee par Texpressidn 

suivante: • N * 

"... <c(t 1 _i)>+g(t 1 )> 

AP m ax(Ti+l)== APmac(Ti). 2.<qii)> ' :. ' 

ou <C(Tj_i)> etWe(Tl)> representent respectivemerit la valeur moyenne 

' du: coOtl«pbur);yfe 
';i5^\'pfeced^t^ 

dans un voisinage born6 du p^ premier temps 

tous les parametres simultan§ment puis en modifiant successiyement un 

• ' parametre a la fojs; .:;::vv-:(-^ : -'-'- • • • .'•:/■• ^.f^/,{- ; ij ; ^^;: r -'^;;^- : -i- •■. ■ 
20 8. Dispdsitiff pb|jr ila mise eh oeuvre du prbc^dd selpn les 

revendications 1 a 6; c^ract§nse en ce qu'il c»mporte des moyens de 
traitement (25) pour la detection tfobjets dans un lot <f images determine, 
parametres par un jeu determine de parametres (P), des moyens de calcul 
(26) de la fonction de cout (C) recevant sur une premiere entree les mesures 

25 delivrees par les moyens de traitement (25) et sur. une deuxieme entree les 
mesures issues de^la yentMerrain, et des mby^is #);<le minimisation du 
■■ cout '(C) recevant; des^er^rt^^ \ sortie un 

nouveau jeu de parametres de reglage des moyens de traitement (25). 

9. Application du procede selon I'une quelconque des 

30 revendications 1 a 7, a la detection d'bbjets ponctuels, caracterise en ce que 
les vecteurs de mesure representent les coprdonnees des objets deletes et 
en ce que I'eiTeur maximale autcfis^ de grandeur de la 

resolution d'un capteur determine^ -'utilise dans le traitement de detection, et; 
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en ce que seuls les taux de fausses alarmes (P fa ) et de non-detections 
(P nd ) sont pris en compte dans le calcul du cout (C). 

10. Application selon la revendication 9, caracterise en ce que le 
traitement comporte des moyens de filtrage pour la reduction du bruit, des 
moyens de seuillage de luminance, des moyens de recherche des 
composantes connexes, des moyens de calcul de la taille des objets 
detectes, et des moyens de selection d'objets dont la taille est comprise 
entre deux valeurs determinees. 
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